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Abstract 

We report a measurement of the charge asymmetry for same-sign dileptons in B°-B° mixing, 
A s \. The data were collected with the Belle detector at KEKB. Using a data sample of 78 fb _1 
recorded at the T(45) resonance and 9 fb _1 recorded at an energy 60 MeV below the resonance, 
we measure A si = (-1.1 ± 7.9(stat) ± 7.0(sys)) x 10" 3 . 

PACS numbers: 11.30.Er, 12.15.Ff, 13.20.He, 14.40.Nd 



*on leave from Nova Gorica Polytechnic, Nova Gorica 
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I. INTRODUCTION 



The Standard Model allows CP violation in B°-B° mixing [1]. In particular, there is a 
possible difference between the B° — > B° and B° — > B° transition rates that can manifest 
itself as a charge asymmetry in same-sign dilepton events in T(AS) decays when prompt 
leptons from semileptonic decays of neutral B mesons are selected. With the assumption of 
CPT invariance, the flavor and mass eigenstates of the neutral B mesons are related by 

\B H )=p\B°) + q \B°), 

\B L }=p\B°}-q\B }, (1) 



where \p\ + \q\ = 1. The time- dependent decay rate for same-sign dileptons is given by 



r T (45H« + (At) = ^e-l A 'l/^ I [cosh At) 



(2) 



for the £ + £ + sample. For the £~£" sample, p/q is replaced by q/p. Here Am^ and Ar are 
the differences in mass and decay width between the two mass eigenstates, t b o is the average 
lifetime of the two mass eigenstates, and At is the proper time difference between the two 
B meson decays. In this analysis only the absolute value of At is used. It is assumed that 
the semileptonic decay of the neutral B meson is flavor specific and CP conserving, so that 
A e = A e , where A e = {X~£ + v £ \B°) and A e = (X + £-u e \B°). If CP is not conserved in 
mixing, the condition \p/q\ = 1 is no longer true and the decay rates for £ + £ + and £~£~ 
samples can differ. As can be seen in Eq. 2, the At dependence is the same for the £ + £ + 
and £~£~ samples, and therefore the CP violation shows up as a At-independent charge 
asymmetry, defined as 

A _ r T (4SHi+i+ - r T (4SHi-i- _ 1 - \q/p\ 4 ^ 4Re(e B ) , . 

Sl r T(4S) ^+ £ + + r T (4S)_^-/- 1 + \q/p\ 4 ~ 1 + |e B | 2 ' 

Here e# corresponds to the ex describing CP violation in the neutral K meson system. 
Standard Model calculations give the size of this asymmetry to be of the order of 10~ 3 [2, 3]. 
A significantly larger value would therefore be an indication of new physics. 

Experimentally, a measurement of same-sign dilepton events that originate from B°B° 
and B°B° initial states requires careful charge-dependent corrections, which are done in 
several steps. First, the contribution from continuum e + e~ — > qq (where q = u, d, s or c) 
to same-sign dilepton events is subtracted using off-resonance data. Second, all detected 
lepton tracks are corrected for charge asymmetries in the efficiencies for track finding and 
lepton identification, and for the probabilities of misidentifying hadrons as leptons. After 
these corrections, the remaining same-sign dilepton events still contain backgrounds from 
B°B° and B + B~ events. The last step of this analysis is to separate the signal events from 
these background events using their different behavior in the At distributions. 



II. BELLE DETECTOR 



This analysis is based on a data sample of 78 fb 1 at the T(4S') resonance ("on- 
resonance") and 9 fb" 1 at 60 MeV below the T(4S') resonance ("off-resonance") collected 
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with the Belle detector [4] at the KEKB asymmetric e + e _ collider [5]. The Belle detec- 
tor is a large-solid-angle magnetic spectrometer that consists of a three-layer silicon vertex 
detector (SVD), a 50-layer central drift chamber (CDC) for tracking, a mosaic of aerogel 
threshold Cherenkov counters, time-of-flight scintillation counters (TOF), and an array of 
CsI(Tl) crystals for electromagnetic calorimetry (ECL) located inside of a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect K\ mesons and to identify muons (KLM). 

A. Track Finding Efficiency 

The track finding efficiency is determined by analyzing a sample where simulated single 
electron or muon tracks are overlaid on hadronic events taken from experimental data. 
Lepton tracks are generated to cover the region of 1.2 GeV/c < p* < 2.3 GeV/c and 30° < 
#i a b < 135°, where p* is the lepton momentum in the e + e~ center-of-mass (cm.) frame and 
#i a b is the angle of the lepton track with respect to the z-axis in the laboratory frame. The z- 
axis passes through the nominal interaction point, and is anti-parallel to the positron beam 
direction. The track finding efficiencies are obtained for each bin of p* (1.2-1.3, 1.3-1.4, 
1.4-1.5, 1.5-1.6, 1.6-1.8, 1.8-2.0 and 2.0-2.3 GeV/c), and d Uh (30°-37°, 37°-50°, 50°-77°, 
77°-82°, 82°-lll°, 111 -119°, 119°-128° and 128°-135°). Figure 1 shows fl lab -averaged track 
finding efficiencies for positive and negative tracks separately and their fractional differences 
as functions of p* for electron and muon tracks. Events in all 9^ regions are combined in 
these plots. The charge dependence of the track finding efficiency for both electrons and 
muons is less than 1.0%. 

B. Lepton Identification 

The most important contribution to the electron identification comes from examination 
of the ratio of the ECL cluster energy to the track momentum measured in the CDC. This 
information is combined with the shower measurement in the ECL, the specific ionization 
measurements (dE/dx) in the CDC, matching information between the ECL cluster position 
and extrapolated position of the CDC track, and the ACC light yield to form an electron 
likelihood C e [6]. 

The two-photon process e + e~ — > (e + e~)e + e - is used to estimate the electron identification 
efficiency. For this data sample, events are required to have: i) two tracks with particle 
identification information inconsistent with a muon hypothesis, laboratory momenta greater 
than 0.5 GeV/c and transverse momenta greater than 0.25 GeV/c; ii) at least one ECL 
cluster with energy greater than 20 MeV. The two tracks are required to have: i) an 
acolinearity angle whose cosine is greater than —0.997; ii) a transverse-momentum sum less 
than 0.2 GeV/c; iii) a longitudinal momentum sum of less than 2.5 GeV/c in the cm. frame; 
iv) an invariant mass less than 5 GeV/c 2 . In addition, the sum of the ECL cluster energies 
has to be between 0.6 GeV and 6.0 GeV. The electron identification efficiency is obtained 
by taking the ratio of the number of tracks selected with the above requirements with and 
without additional electron identification requirements. 

For muon identification, CDC tracks are extrapolated to the KLM and the measured 
range and transverse deviation in the KLM is compared with the expected values to form 
a muon likelihood £ M [7]. The muon identification efficiency is determined by analyzing a 
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FIG. 1: #i a b-averated track finding efficiencies as a function of cm. momentum for positron 
tracks e+ (a), electron tracks e_ (b), and charge dependence denned as (e+ — e_)/(e+ + e_) (c). 
Corresponding plots for muon tracks are shown in (d), (e), and (f) 

data sample where simulated single-muon tracks are overlaid on the hadronic events taken 
from experimental data. 

The lepton identification efficiencies are obtained in the same p* and #i a b bins used for 
the track finding efficiency study. Figure 2 shows the 6*i a b-averaged charge-dependent lep- 
ton identification efficiencies, where electron tracks are required to satisfy C e > 0.8, and 
the muon tracks are required to satisfy > 0.9 and to have a reduced x 2 value for their 
transverse deviations in the KLM of less than 3.5. The charge dependence of the identifica- 
tion efficiency is less than 1% for both electrons and muons. Exactly the same criteria are 
imposed to select signal leptons. 

C. Hadron Misidentification 

The information on the ACC light yield, the time-of-flight measurement, and dE/dx 
measurements in the CDC is combined to provide hadron likelihoods: C w for pions, Ck for 
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FIG. 2: The #i a b-averaged identification efficiencies as a function of cm. momentum for positrons 
e+ (a), electrons e_ (b), and charge dependence defined as (e+ — e_)/(e+ +£-) (c). Corresponding 
plots for muons are shown in (d), (e), and (f). 

kaons, and C p for protons. 

The hadron fake rate, which is defined as the probability that a hadron track is mis- 
identified as a lepton, is determined from a sample of K s — > 7r + 7r~ decays for pions, 
— > K + K~ decays for kaons, and A — > £>7r~ (A — > j>7r + ) decays for protons. These de- 
cays are selected from a hadronic event sample described below. To select these track-pair 
combinations, the distance of closest approach with respect to the run-dependent interaction 
point and the position of the decay vertex are used. The difference in z position of the two 
tracks and the angle between position vector of the decay vertex and the momentum vector 
of K s or A(A) in r-0 plane at the decay vertex are also used. When evaluating the fake 
rate for positively charged tracks, invariant masses are plotted for a sample of pairs where 
both tracks satisfy the above criteria and the negative track satisfies an additional hadron 
identification requirement. The signal yields are obtained by fitting the resulting mass dis- 
tribution to a sum of a double Gaussian signal and a smooth background function in two 
ways: once without imposing any lepton identification requirement and once after imposing 
the lepton identification requirement on the positive track. The ratio of the two signal yields 
gives the fake rate for the positive charged tracks. The following requirements are imposed 
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on the likelihood ratios: C n / (C n + Ck) > 0.8 for pions in Kg — > 7r + 7r~, Ck/{C>k + £?r) > 0.8 
for kaons in — > K + K~, and C n /(C n + £ p ) > 0.8 for pions in A — > pn~ . The fake rate 
for negatively charged tracks is obtained by repeating the procedure above with the roles of 
negative and positive tracks reversed. 

In this study, because of low statistics, #i a b-averaged values are used, and some p* bins 
are combined. 

Figure 3 shows the hadron fake rates as a function of momentum in the laboratory frame. 
The rate of pions faking electrons is at most 0.1% for both charges. The rate of kaons faking 
electrons decreases rapidly as pi a b becomes larger and is less than 0.2% for pi ab > 1.4 GeV/c, 
with no significant charge dependence. While the rate of protons faking positrons is nearly 
zero, the rate for anti-protons faking electrons is as large as 4% due to the large anti-proton 
annihilation cross section in the ECL. Because of low statistics, the rate of protons faking 
positrons is obtained from the ^b-averaged value over all momenta. The rate of pions faking 
muons is about 1% for p\ a ^ > 1.5 GeV/c and shows no significant charge dependence. The 
rate of kaons faking muons is 1% to 2%, with that for K + being about 50% larger than that 
for K~ due to the larger kaon-nucleon cross section for the K~ . The rate of protons faking 
muons is less than 0.4% and shows no clear charge dependence. 

III. EVENT SELECTION 

A. Hadronic Event Selection 

Hadronic events are required to have at least five tracks, an event vertex with radial 
and z coordinates within 1.5 cm and 3.5 cm, respectively, of the nominal beam interaction 
point, a total reconstructed cm. energy greater than 0.5 W (W is the cm. energy), a 
net reconstructed cm. momentum with a z component less than 0.3 W/c, a total energy 
deposited in the ECL between 0.025 and 0.9 W, and a ratio R2 of the second and zeroth 
Fox- Wolfram moments [8] of less than 0.7. 

B. Dilepton Event Selection 

Lepton candidates are selected from among the charged tracks by requiring the criteria 
previously described. In both electron and muon cases, a distance of closest approach to the 
run-dependent interaction point less than 0.05 cm radially and 2.0 cm in z is required. At 
least one SVD hit per track in the r-<f> view and two SVD hits in the r-z view are required. 
To eliminate electrons from 7 — > e + e~ conversions, electron candidates are paired with all 
other oppositely charged tracks and the invariant mass (assuming the electron mass) M e + e - 
is calculated. If M e + e - < 100 MeV/c 2 , the electron candidate is rejected. If a hadronic event 
contains more than two lepton candidates, the two with the highest cm. momenta are used. 

The two lepton candidates must satisfy additional criteria. The cm. momentum of each 
lepton is required to be in the range 1.2 GeV/c < p* < 2.3 GeV/c. The lower threshold 
reduces contributions from secondary charm decay; the upper threshold reduces continuum 
contributions. Each lepton track must be in the range 30° < 9^ < 135°, where it has better z 
vertex resolution and lepton identification. Events that contain one or more J/ ip candidates 
are rejected. The invariant mass of each candidate lepton paired with each oppositely 
charged track (assuming the corresponding lepton mass) is calculated. If the invariant mass 
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FIG. 3: Rates of pions ((a) and (d)), kaons ((b) and (e)), and protons ((c) and (f)) faking electrons 
and muons vs. laboratory momentum. Filled circles are for positive tracks and open circles are for 
negative tracks. The increase in the rate of kaons faking electrons at low momentum clearly visible 
in (b) is due to the overlap of the electron and kaon energy-loss bands; the other distributions are 
discussed in the text. Note that the minimum values of vartical axes are set to negative values. 
The dotted lines shows fake rate equal to zero. 

falls within the J/ip region, denned as -0.15 GeV/c 2 < (M e + e - - Mj^) < 0.05 GeV/c 2 or 
-0.05 GeV/c 2 < (M M+M - -Mj/^) < 0.05 GeV/c 2 , the candidate event is rejected. The looser 
lower mass window for the electron pair invariant mass accomodates bremsstrahlung of the 
daughter electron(s). 

As can be seen in Figure 4, distributions of the opening angle of the two tracks in the 
cm. frame, cos 9* u , for the \x\x and e\x pairs show distinct peaks in the back-to-back direction 
(cos#^ ~ — 1). This background is caused by hadron tracks misidentified as muons among 
jet-like continuum events. Also, spikes can be seen at cos6* u = 1 for the fxfj, pair events. This 
structure is caused by jet-like continuum events where a non-muon track is identified as a 
muon because hits in the KLM from a nearby true muon are assigned to it. The dilepton 
opening angle in the cm. frame 6* u is required to satisfy —0.80 < cos^ < 0.95 in order to 
reduce this background. 

With these selection criteria there are 46533 positive and 45477 negative same sign 
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FIG. 4: cos#^ distributions for dilepton samples in the on-resonance (open histogram) and scaled 
off-resonance (filled histogram) data. (a),(b) show ee events, (c),(d) show [X[i events and (e),(f) 

are from e// combinations. (a),(c) and (e) are the ++ charge case and (b), (d) and (f) are the 

charge case. The arrows indicate the selection requirements. 



dilepton events found in the on-resonance data. Continuum contributions are estimated 
to be 2230 for positive and 1574 for negative same sign events, based on the yield from 
off-resonance data. To estimate the continuum contribution from off-resonance data, off- 
resonance yields were scaled by the integrated luminosities and cross-section ratio. The 
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scaling factor is defined by 

/= J^** (4) 

where J C on ^ g)dt and s n(off) are the integrated luminosities and the square of cm. energies 
for on(off)-resonance, respectively. 

These dilepton yields, decomposed into three lepton categories, are given in Table I. 



TABLE I: Numbers of dilepton events. 





On-resonance 


Off-resonance 


Continuum 


Combination 


positive 


negative 


positive 


negative 


positive 


negative 


ee 


9059 


9028 


11 


11 


96.2± 28.9 


96. 2± 28.9 


nn 


14672 


14014 


144 


100 


1259.2±104.9 


874.4± 87.4 




22802 


22435 


100 


69 


874.4± 87.4 


603.4± 72.6 


total 


46533 


45477 


255 


180 


2229.8±139.6 


1574.0±117.3 



C. Az Determination 

The ^-coordinate of each B meson decay vertex is the production point of the daughter 
lepton, which is determined from the intersection of the lepton track with the run-dependent 
profile of the interaction point. The distance between the ^-coordinates of the two leptons, 
\Az\, is defined as \Az\ = \z(£i) — z(£ 2 )\- 

In order to estimate the effect of detector resolution in the Az distribution, J/ip decays 
to e + e~ and fi + fi" are used. In these events, the two tracks originate from the same point, 
so the measured Az, after the background contribution is subtracted, yields the detector 
resolution. Candidate J/ ip mesons are selected using the same requirements used for dilepton 
events, except for the J/ip veto. The J/ip signal regions are defined as 3.00 GeV/c 2 < 
M(e + e") < 3.14 GeV/c 2 and 3.05 GeV/c 2 < M(n+n~) < 3.14 GeV/c 2 and the sideband 
region as 3.18 GeV/c 2 < M(£ + £~) < 3.50 GeV/c 2 for both electrons and muons. 

The invariant mass distributions of J / ip candidates are fitted to a function given by 

(M-M f (M-Mp) 2 

N(M) = h e + h x e ^ + A(M - B) 2 + C. (5) 

Here, h and hi are the normalizations of the two Gaussians used to describe the signal, 
M is the Gaussian mean which is common to both Gaussians, and a and 04 are the 
standard deviations of the Gaussians. A parameter S, defined as S = cr for M > M and 
S = a + ct(M — M ) for M < M , is introduced to modify the lower mass tail of one of the 
Gaussians to take the effect of bremsstrahlung into account. A, B and C are the parameters 
of the background function. The Az distribution of the sideband region is scaled to the 
background yield in the signal region and subtracted from the signal region Az distribution. 

The J/ip mass distributions and the Az distributions are shown in Figure 5. The RMS 
values of the Az distributions are 193 pm for J/ip — > e + e", 177 pm for J/ip — > and 
185 pm for the combined J/ip — > £ + £~ sample. 
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FIG. 5: Mass distributions for J/ip — > e + e~ (a) and J/^ — > (b). The arrows indicate the 

selection criteria for each decay mode. The dashed lines indicate the fitted background compo- 
nent and the solid lines show the total fit results. The cross hatched area shows the estimated 
combinatorial background in the signal region and the single hatched area is the sideband region 
used to estimate the background Az distribution. The Az distributions for J/ijj — > e + e~ (c) and 
J/tjj — > fJ, + fx~ (d). Open histograms are for all J/tp candidates in the signal region and hatched 
histograms are for the background. 
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D. Subtraction of Continuum Events 

A sample of dilepton events originating from BB events is obtained by subtracting the 
luminosity and cross-section scaled off-resonance data from the on-resonance data. Since 
the kinematics of dilepton candidates in BB decays is different from those in continuum 
events in each of the variables (p*, p%, 0%, 0* u , Az), where 6\, 2 s is the polar angle in the 
cm. frame with respect to the beam axis of the more (less) energetic lepton, the subtraction 
should, in principle, be performed in this six-dimensional space, separately for each lepton 
flavor and charge combination. 

Given the available statistics, this approach is not possible. Instead, we perform the sub- 
traction by weighting the on-resonance and off-resonance yields for one of the six kinematical 
variables, while integrating over the five other variables. We obtain weighting factors for the 
six variables by repeating this procedure. Since, to a first approximation, the six variables 
are not correlated with each other, this approach provides the BB yield in the six variable 
space. The weighting factors are given by w(k) = {l/rBB){N on (k) — fN g(k))/N on (k) where 
k denotes each of six variables, / is the scaling factor for the luminosity and cm. energy 
introduced in Eq. 4, and vbb = ATj^ al /iv"*° tal is the fraction of total BB events in the on- 
resonance yield after integrating over all six variables; it is used for normalization. While 
the weighting factors show very little dependence on p\, p 2 , 0*, and 0\ for all combinations 
of lepton flavors and charges, a clear dependence is observed for 0* u in the case of the \x\x 
and e/i data samples as shown in Figure 6. A clear dependence on Az is also seen for all 
lepton pair combinations. 

Using this method, the dilepton candidate yield for each lepton flavor and charge combi- 
nation is then given in terms of the on-resonance yield and the weighting factors by 

N BB (P1,P*2, 0{, 9* 2 , 0h A*) = r BB Y[w(k)N on (pl,P*2, 0*i, ^ (6) 

k 

The Az distributions of the dilepton yields are obtained by projecting 
Nbb(.Pi, P2, &2, &u, Az) onto the \Az\ axis. 



IV. RESULT 

A. Corrections to Lepton Candidates 

The number of detected leptons for each lepton flavor and charge N^ et is related to the 
number of true leptons Nf by 

JV£t(p*, 01ab) = Nftp\ 01ab)4k(P*, M | S%(p\ Iab ) + r ^*> *tabto«(P*, ftab) } ,(7) 

I h=ir,K,p ) 

where and 5p id are the efficiencies for track finding and lepton identification, is the 
relative multiplicity of hadron of type h with respect to leptons of type £ in BB event, and 
77^ is the rate of hadrons h faking leptons I. The relative multiplicities are determined from 
BB Monte Carlo (MC) events, and are shown in Figure 7. The relative multiplicities include 
the effects of decay-in-flight and interaction with the detector. 
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FIG. 6: cos 6* u dependence of the weighting factor for the fraction of the dilepton yield originating 
from BB decays in the on-resonance data for e + e + (a) and e~e~ (b) and corresponding quantities 
for nn (c) and (d) and e/j, (e) and (f). 

Using the measured efficiencies, fake rates and the MC-determined relative multiplicities, 
the correction factors Ng/N&n, are determined in 7 bins of p* (1.2-1.3, 1.3-1.4, 1.4-1.5, 
1.5-1.6, 1.6-1.8, 1.8-2.0 and 2.0-2.3 GeV/c), and 8 bins of 9 lab (30°-37°, 37°-50°, 50°-77°, 
77°-82°, 82°-lll°, m°-119°, 119°-128° and 128°-135°). (The fake rates are measured in 
the laboratory frame, but are converted into p*- and ^independent measurements for this 
correction.) 

After the correction, the dilepton sample contains true leptons that come either from 
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FIG. 7: Relative multiplicities of hadrons as a function of pi a b with respect to leptons. (a), (b) and 
(c) are for electrons and (d), (e) and (f) are for muons. (a) and (d) are pion to lepton, (b) and (e) 
are kaon to lepton, (c) and (f) are proton to lepton ratios. Filled circles are for positive tracks and 
open circles are for negative tracks. The effects of decay-in-flight and interaction with detector are 
included. Though in (a), (b), (d) and (e), the difference between positive and negative cases is less 
than 1%, in both (c) and (f), the proton rate is larger than the anti-proton rate by about 8%. The 
increases in multiplicities for pions and kaons around 3.5 GeV/c are due to two-body hadronic B 
decays. 

prompt neutral B meson decay (signal) or from background processes such as charged B 
meson decay, secondary charm decay, or other leptonic B meson processes. 



B. Fit to Az Distribution 

A binned maximum likelihood fit with signal and background contributions is used to 
extract A s i(|A,2:|) from the Az distribution. The overall background level is obtained from a 

fit with positive (++) and negative ( ) samples combined. In this fit the Az distribution 

for signal events is given by Eq. 2, assuming Ar is small [9], 

p ss K e -|At|/r fl0 _ cos (Am ( iAt)), (8) 
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convolved with the detector response function described earlier. Here, t b o and Anid are 
fixed to their world average values [9]. 

The backgrounds are divided into two categories: correctly tagged (CT), and wrongly 
tagged (WT). The CT category mainly contains events in which both leptons come from 
secondary charm decay in B°B° — > B B°(B°B ) (mixed) processes. The WT category con- 
tains events in which one lepton is from secondary charm decay of unmixed B°B° or B + B~ 
and the other is from a semileptonic B decay. Though background Az distributions are 
estimated using MC simulations, the MC underestimates the width of the Az distribution. 
To correct for this, the MC Az distribution is convolved with a Gaussian of standard devi- 
ation a = 69 [im [10]. The Az distribution for the true same-sign dilepton events, where 

positive (++) and negative ( ) samples are combined, is shown in Figure 8- (a) together 

with the fit results. The x 2 /n.d.f '. of the fit is 48.75/38 (note that only statistical errors are 
included). In the fit, the ratio of CT to WT is fixed at the MC value, and only the ratio of 
signal and background is allowed to float. The MC-estimated CT and WT contributions to 
the Az distribution are shown in Figure 8-(a). 



C. Charge Asymmetry 

The measured same-sign dilepton charge asymmetry is defined as 

A (\A N++{Az)-N-{Az) 

MAZ) = N++(Az)+N-(Azy (9) 

where iV ±± (Az) are the true dilepton yields as a function of Az. 

Since N ±:t (Az) are the sum of signal and background, N ±:iz (Az) = (Az) + N b ±i: (Az) , 
the dilepton charge asymmetry A s \ is related to Am by 

, , N++(Az) - N--(Az) NJAz) . .. . I/A . 

Aee(Az) = s v ' . - .* v — - /1tr/A x J /A xx = AJAz)d(Az) 10 
v ; N s (Az) (N s (Az) + N b (Az)) K J K ' K ' 

where N s = N+ + + N~~ and N b = + N b ~~ . Here, N b ++ = N b " is assumed. A 

dilution factor, d(Az) = N s (Az) / (N s (Az) + N b (Az)), is calculated using the signal and 
background yields, which are determined in the fit given in Figure 8-(a). The result for 
A s i(Az), determined from the measured A u (Az) and d(Az), is shown in Figure 9. 

The dilepton charge asymmetry A s \ is a time integrated quantity and does not depend 
on Az. Fitting this distribution to a constant in the region of 0.015 cm < \Az\ < 0.200 cm, 
yields A s \ = (— 1.1±7.9) x 10~ 3 , with a x 2 / n -d-f- of 36.20/36 (again note that only statistical 
errors are included). The optimum fitting range is determined using a MC study that 
maximizes N s /y/N s + N b . 



D. Cross Checks 

As a consistency check, A s \ is obtained separately for the ee, /x/x, and e/x data samples. 
The results, A sl (ee) = (-13.3 ± 14.4) x 10" 3 , A sl (fifi) = (+16.5 ± 17.3) x 10~ 3 , and A sl (efx) = 
(—3.6 ± 11.3) x 10 -3 , are consistent with the primary result. Here the errors are statistical 
only. 
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FIG. 8: (a) Az distribution for the true dilepton events (++ and are combined). Points with 

error bars are data. The dot-dashed line shows the contribution from CT backgrounds, the dotted 
line shows the WT background contributions, the dashed line indicates the signal component and 
the solid line indicates the fitted total, (b) Difference between data and fit result as a function of 
|Az|. 



The validity of the assumption, to extract A s \ from Au(Az), is confirmed by repeating 
the fit without it. This yields iV+ + = 20200 ± 212 and N~~ = 19766 ± 210 in the range 
0.015 cm < \Az\ < 0.200 cm, which is consistent with the initial assumption. 



17 



■a 0.5 



0.4 



0.3 



0.2 



0.1 



i 



r v + ' 



-0.1 



-0.2 



-0.3 



-0.4 



-0.5 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

lAzl(cm) 



FIG. 9: |Az| distribution for A sl . 



E. Systematic Errors 



Systematic errors in the determination of A s \ come from uncertainties in: i) the event se- 
lection criteria; ii) the continuum subtraction; iii) corrections for efficiencies of track finding 
and lepton identification and for lepton misidentification; iv) the dilution factor determina- 
tion; and v) the Az fit for the determination of A s \. 

Uncertainties in the selection criteria are estimated by repeating the analysis with differ- 
ent threshold values. For the track selection, the #i a b lower (higher) limit is varied from the 
nominal 30° (135°) to 40° (125°) in ten (eighteen) steps, the closest-approach criterion in r-0 
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from its nominal value of 0.05 cm to 0.02 cm in fourteen steps, the closest-approach criterion 
in r-z from 2 cm to 1 cm in fourteen steps, the p* lower (higher) limit from its nominal value 
of 1.2 GeV/c (2.3 GeV/c) to 1.3 GeV/c (2.2 GeV/c) in fourteen steps, the requirement on 
the number of SVD hits by ±1 from its nominal value of greater than or equal to two for 
r-z and one in r-<p. For the event selection, the cos 9* u lower (higher) criterion is varied from 
nominal —0.80 (+0.95) to —0.60 (+0.70) in twelve (seventeen) steps. In addition, the mass 
windows that reject J /if) and 7 — > e + e~ are widened by 20% and 50%, respectively. 

To estimate the systematic error from the continuum subtraction, the analysis is repeated, 
varying the off- resonance yield by ±ler. 

For systematic uncertainties arising from the corrections to lepton candidates, we consider 
separately the effects due to the uncertainties in each bin and possible correlated sources of 
error. The fake rates from pions, kaons and protons and the lepton identification efficiencies 
are assumed to be independent. The systematic uncertainties from these track corrections 
are estimated by varying the efficiencies and fake rates by ±1(7 bin by bin. The system- 
atic effects associated with the relative hadron multiplicities are estimated by varying the 
values by ±ler for all bins (both electron and muon) simultaneously. The uncertainties for 
each hadron species are estimated separately, assuming that the multiplicities of each are 
uncorrelated. 

The contribution from the detector Az response function is estimated by changing the 
response function shape according to the statistics of each Az bin of J/ip — > £ + £~ sam- 
ple. The contribution from the 69 [im smearing is estimated by repeating the analysis with 
43 [im and 96 [xm smearing, corresponding to the uncertainty in the difference of the res- 
olution between data and MC. The contributions from uncertainties in Am^ and r B o are 
also estimated by varying the nominal values by ±ler. The fitting range of the dilution 
factor determination is varied from its nominal range of 0.000 cm < \Az\ < 0.200 cm to 
0.025 cm < \ Az\ < 0.200 cm. For the fitting range for the determination of the final A 8 \, the 
lower limit is varied from 0.000 cm to 0.045 cm in nine steps. The results of the systematic 
error determination for A s \ are summarized in Table II. 



The measured A s \ is consistent with zero, or equivalently, \q/p\ is consistent with unity. This 
implies that CP violation in B°-B° mixing is below the O(10 -2 ) level. The CP violation 
parameter e# can be calculated as 



V. CONCLUSION 




Re(e B ) 

1 + M 2 



(-0.3 ± 2.0(stat) ± 1.7(sys)) x 10~ 3 , 



(13) 



using the (exact) formula 



Rejeg) 

1 + M 




1 - y/(l - A sl )/(1 + A al ) 

1 + v/(i-A sl )/(i + A sl )' 



(14) 
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TABLE II: Source of systematic errors for the measurement of A s \. 



Category 


Source 


A A s i (xl(T 3 ) 


Event selection 


Track selection 


±2.61 




cos9* u cut 


±0.63 




Lepton pair veto 


±2.33 


Continuum subtraction 




±4.88 


Track corrections 


Track finding efficiency 


±1.56 




Electron identification efficiency 


±0.56 




Muon identification efficiency 


±1.98 




Fake electrons 


±0.45 




Fake muons 


±0.81 




Relative multiplicity 


±0.56 


Az fit for dileptons 


Detector response function 


±0.07 




Am d 


±0.08 






±0.07 




69 um smearing of background Az 


±0.13 




Statistics of signal MC 


±0.01 




Statistics of background MC 


±0.19 




Fitting range 


±0.04 




Assuming N^~ + = N b 


±1.59 


Az fit for A sl 


Fitting range 


±1.30 


Total 




±6.97 



These results are consistent with previous measurements [11] and provide significantly more 
restrictive bounds. Previous measurements of Re(ee)/(1 + |eg| 2 ) are listed in Table III 
together with this measurement. 

TABLE III: Comparison of the Re(e.g)/(1 + I^bI 2 ) measurements. 



Experiment 


Re(e B )/(l±|e B | 2 ) 


ALEPH 


-0.003 ± 0.007 


CLEO 


0.0035 ± 0.0103 ± 0.0015 


BABAR 


(1.2±2.9±3.6)xl0" 3 


This experiment 


(-0.3 ± 2.0 ± 1.7) xlO -3 
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